Entanglement in a Raman-driven cascaded system 
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The dynamics of a cascaded system that consists of two atom-cavity subsystems is studied by 
using the quantum trajectory method. Considering the two atom-cavity subsystems driven by a 
Raman interaction, analytical solutions are obtained. Subsequently, the entanglement evolution 
between the two atoms is studied, and it is shown that the entanglement can be stored by switching 
off the Raman coupling. By monitoring the radiation field, the entanglement between the two atoms 
can be enhanced. 
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I. INTRODUCTION 

The concept of entanglement has been of great inter- 
est since the early days of quantum mechanics [1|, and 
it has become of central importance in a variety of dis- 
cussions on the fundamental aspects of the theory [3, Q ■ 
Nowadays entanglement is receiving new attention in the 
context of the rapidly developing fields of quantum infor- 
mation, quantum computation and quantum technology; 
for reviews, see [j, [a, |a, 0] ■ Entanglement is considered 
to be the characteristic feature that allows quantum in- 
formation to overcome some of the limitations imposed 
by classical information. Cold trapped atoms interact- 
ing with quantized light fields are promising candidates 
for the realization of quantum computing and quantum 
communication protocols 8, 3] . The combination of long- 
lived atomic states and light fields can be used in quan- 
tum networking for the distribution and processing of 
quantum information [lO . Ill| . In the context of entan- 
glement preparation between atoms at separate nodes, 
a variety of schemes have been proposed, for example, 
by measuring the superpositions of light fields released 
from separate atomic samples, or by measuring a probe 
light field that has interacted in a prescribed way with 
different samples. Due to the indistinguishability in the 
measurement, and conditioned on the results of the mea- 
surements, the atomic system is projected into an entan- 
gled state [l2, [l3, llJ, lUl • An unconditional preparation 
of entanglement has also been analyzed in the case of 
a cascaded system. This unconditional preparation has 
been discussed for two distantly separated atoms pia.lTn|. 
as well as for separate atomic ensembles [l^. Moreover, 
the recent achievements in cavity QED and in tapped ion 
techniques have rendered it possible to experimentally 
generate pairs of entangled atoms [l9|, to create entan- 
gled states of several atoms [20| , and even long-lived en- 
tanglement of two macroscopic ensembles of atoms [2l| . 

In the spirit of these previous achievements, in the 
present contribution we will consider the quantum trajec- 
tory approach for a cascaded open quantum system [22| . 
We study the dynamics of a system that consists of two 
atom-cavity sub-systems A and B. The quantum source 
A emits a photon and the second quantum subsystem B 
reacts on the emitted photon. We will first consider an 



unconditional preparation of the entanglement between 
the two atoms. Second, the effects of a null-measurement 
conditional preparation is analyzed with respect to a pho- 
todetector of a given efficiency monitoring the field radi- 
ated by the cascaded system. 

As is clearly discussed in Refs. [la, [ijjilal' ^^^ advan- 
tage in using a cascaded system is that the dynamical 
evolution of the open quantum system itself creates the 
entanglement. It is an unconditional preparation, and 
it is not related to a "click" or "no click" at a detector, 
where the measurement projects the atomic system onto 
the desired entangled state. In this sense, we can say 
that it is a dynamical generation of entanglement, and 
not a conditional one. One could also use a detector of 
given efficiency to monitor the radiated field, to prepare 
the system conditioned upon "no click" at the detector. 
This allows us to realize a quantum state preparation 
conditioned upon the limited knowledge of the observer's 
imperfect detector. In this way one can combine the ad- 
vantages of using a cascaded system, with its intrinsic 
dynamical generation of entanglement, and a conditional 
preparation with a detector of non-unit efficiency. In 
the case under study the entanglement between the two 
atoms can only increase due to this conditional prepara- 
tion, even for imperfect detection. In the limiting case 
of a detector of zero efficiency, we return to the case of 
unconditional preparation, where only dynamically gen- 
erated entanglemet is present. In the system under study 
we consider a Raman configuration for driving the atom- 
cavity interaction. By switching off the lasers beams, 
the Raman interaction vanishes, so that the entangle- 
ment generated between the atoms remains unchanged 
and can be stored. 

The paper is organized as follows. In Sec.[n]the master 
equation describing the dynamics of the cascaded system 
is introduced, and the problem is solved analytically by 
using the quantum trajectory method. In Sec. IIIII the 
unconditional preparation and storage of the entangle- 
ment between the two atoms is analyzed. The condi- 
tional preparation of the entanglement between the two 
atoms is discussed in Sec. |TVl Finally, some concluding 
remarks are given in Sec. |Vl 



II. CASCADED SYSTEM DYNAMICS 

In this section we analyze the dynamics of the system 
under study. The cascaded open quantum system con- 
sists of two atom-cavity subsystems A and B, where the 
source subsystem A is cascaded with the target subsys- 
tem B, as sketched in Fig. [T] The cavities have three 
perfectly reflecting mirrors and one mirror with trans- 
mission coefficient T ^ 1. In the two subsystems A and 
B, denoted by A; = a, 6, respectively, we consider a three- 
level atom coupled to a cavity mode of frequency tuk via 
a Raman interaction, as indicated in Fig. [21 This config- 
uration is obtained by irradiating the atom with a laser 
beam of frequency w^ such that Wfc — w^ = wfg, where 
uj^Q is the transition frequency between the two atomic 
energy eigenstates |lfe) and |Ofe). The laser beam is de- 
tuned from the electric dipole transition \lk) <-> \2k) by 
Afc, chosen to enhance the Raman-coupling strength, but 
also to avoid electronic excitations. The Rabi frequency 
of the laser is denoted by fife and gk is the strength of 
coupling between the cavity mode and the |Ofe) <-!■ |2fc) 
transition. The cavity mode is damped by losses through 
the partially transmitting cavity mirror. In addition to 
the wanted outcoupling of the field, the atom can spon- 
taneously emit a photon out the side of the cavity, or a 
photon can be absorbed or scattered by the cavity mir- 
rors. 

To describe the dynamics of the system we will use 
a master equation formalism, and solve it by using the 
quantum trajectory method [23, [2J, l2^| ■ For sufficiently 
large detuning, gk/ ^k and fife/Afc <C 1, the excited state 
|2fc) will not become significantly populated and can be 
adiabatically eliminated. This leads, treating the dissipa- 
tion due to the cavity losses in a standard way [26l.[27l.[28j. 
to the following master equation for the reduced density 
operator p(t) of the system: 



dp{t) 



dt ih 
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The Hamiltonian is given by 



H = Ha + Hb+ ih 



\/l^al^b 



-%a) - e'H^ 



(1) 



(2) 



where Ha and Hb describe the atom-cavity interaction 
in the two subsystems A and B, respectively. In the 
rotating-wave approximation they are given by [29| 

Ha = hga (aiio + aUoi) +hl^'^Aii + hKaa) aAoo , (3) 

and 

He = hgb (bBio + b^ Boi) + HA'^Bn + HAbb^bBoo . (4) 



The third term in Eq^J[J|) describes the coupling between 
the two cavities |22l . [30|. In these expressions, d and d^ 
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FIG. 1: The cascaded open system consisting of two Raman- 
driven atom-cavity subsystems A and B. The dashed arrows 
indicate the laser fields needed for the Raman coupling. A 
photodetector PD can be used to monitor the radiation field. 



are annihilation and creation operators for the cavity field 
A, and similarly b and b^ for the cavity field B. We have 



-I] - \ia){ja\ {i,j 



0, 1), and Bi 



Vb){jb 



also defined Ai 
{i,j — 0, 1). In addition, gk = —gk^k/Ak is the effective 
atom-cavity coupling constant and A5, = — ri|/Afe, A^ = 
—g'^/Ak are the two Stark shift terms. Moreover, Ka, Kt 
are the cavity bandwidths and the phase (j> is related to 
the phase change upon reflection from the source output 
mirror, and/or to the retardation of the source due to 
the spatial separation between the source and the target, 

cf. [sii]. 

The jump operators Jt in Eq. ([1]) are defined by 



Ji = V^^a + V^be "''b , (5) 

which describes a photon emission by the cavities; 



J2 



K'aa, 



Js 



k'J), 



(6) 



are associated with a photon absorption or scattering by 
the cavity mirrors; 



and 



Ji = yPal^aAoi+gadAm) 

Jz = ^/T^maAii^^gadAio), 



Jg — V^bip-bBoi+gbbBoQ 
^7 = Jr'MbBn+gbbBw 



(7) 



(8) 



are related to a photon spontaneously emitted by the 
atoms. Here n'f. is the cavities mirrors' absorption and 
scattering rate. Moreover, F/j=7fc/A^ and r'j,=7[,/A^, 
where jk and 7^ are the dipole relaxation rates of the 
atomic state \2k) to the states |Ofc) and |lfc), respectively. 
These relaxation rates are considered to be small in com- 
parison with the detuning. Note that the operator Ji 
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FIG. 2: Raman-type excitation scheme for the atom-cavity 
subsystem A [k = a) and B {k = b). The cavity mode of 
frequency uik and the laser of frequency ojj. are detuned by 
Afc from the atomic state |2fc). The dipole relaxation rates of 
these states to the states |Ofc) and |lfc) are 7fe and 7^, respec- 
tively. 



contains the superposition of the two fields radiated by 
the two cavities, due to the fact that radiated photons 
cannot be associated with photon emission from either A 
or B separately [2a |. 

In the following we will identify, for notational conve- 
nience, the state \a) with the state |1, 0,0,0), which de- 
notes the atom A in the state JIq), the cavity A in the vac- 
uum state, the atom B in the state |0b) , and the cavity 
B in the vacuum state. In the state \h) = |0, 1,0, 0) the 
atom A is in the state |0a) , and the cavity A is in the one- 
photon Fock state. Similarly, we define \c) = |0, 0,1,0), 
\d) = |0, 0,0,1), and |e) = |0, 0,0,0). The state \a) wih 
be considered as the initial state of the system. It fol- 
lows that the Hilbert space that describes the cascaded 
system under study is, in our model, spanned by the five 
state vectors \a), \b), |c), \d), and |e). 

To evaluate the time evolution of the system we use a 
quantum trajectory approach [25, i24i25i|. Note that the 
probability for a jump Ji to occur in the time interval 
[t, t + dt) is given by pi{t) = {J^Ji)t dt. This implies that 
the total probability for a jump due to a spontaneous 
emission in the time interval [t, t + dt) is given, using 
Eqs. 0, ®, ©, and (HI), by 



7 

E 
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{Jjl)tdt 






-1', 



hAh 



{HB)tdt.{9) 



This relation shows that in a time interval ~ A/{gi}) the 
probability to have a jump due to spontaneous emissions 
is ^ (7a+7a+7&+7b)/A, cf. 32]. For a large detuning this 
probability is small. If one is interested to follow the dy- 
namical evolution of the system for several Rabi oscilla- 
tions, in general the effects due to spontaneous emissions 
cannot be neglected [33]. In the present contribution, the 
Raman dynamics would be actually used only for a few 
Rabi oscillations and we may neglect the terms in the 
master equation related to spontaneous emissions. 

Let us now consider the system prepared at time t^ = 
in the state \a). To determine the state vector of the sys- 
tem at a later time t, provided that no jump has occurred 



between time to and t, we have to solve the nonunitary 
Schrodinger equation 



«ftl7lV^no(i)) = H' li'Ut)) 

dt 



(10) 



where H' is the non-Hermitian Hamiltonian given by 



^f'^fU.-ir. ,^^_,ftf^ata 



i 



where we have defined 



ICb — Kb + n'l,. 



(11) 



(12) 



If no jump has occurred between time to and t, the system 
evolves via Eq. (fTO|) into the unnormalized state 



\4'no{t)) = ait)\a) + m\b) + l{t)\c) + S{t)\d) . (13) 

In this case the conditioned density operator for the 
atom-cavity system is given by 



Pnoit) = 



\^no{t)){'^no{t)\ 
{4'no{t)\tpno{t)) 



(14) 



Here we have used the word conditioned to stress the fact 
that this is the density operator at time i, conditioned 
on the fact that no jump has occurred between time to 
and t. 

The evolution governed by the nonunitary Schrodinger 
equation (jlOp is randomly interrupted by one of the three 
kinds of jumps Ji, cf. Eqs. ([5]) and ([6]). If a jump has 
occurred at time ij , ij S {to , t] , the state vector is col- 
lapsed in the state |e) due to the action of one of the 
jump operators. 



J^ l^no(ij)) ^ |e) (^ = l,2,3). 



(15) 



In the problem under study we may have only one jump. 
Once the system collapses into the state |e), the nonuni- 
tary Schrodinger equation ^TU\i lets it remain unchanged. 
In this case the conditioned density operator at time t is 
given by 



Pycsit) = |e)(e| 



(16) 



where we indicate with "yes" the fact that a jump has 
occurred. 

In the quantum trajectory method, the density oper- 
ator p{t) is obtained by performing an ensemble average 
over the different conditioned density operators at time t. 
In the present case, starting at time tg with the density 
operator po — \a){a\, the ensemble average is performed 
over the two possible realizations (histories) "yes" and 
"no" , yielding the statistical mixture 



/5(i) = Pnoit)pno{t) +Pycsit)Pycsit) 



(17) 



Here Pno{t) and Pycs{t) are the probability that between 
the initial time tg and time t no jump and one jump has 
occurred, respectively. Of course, Pno{t) + Pycs{t) — 1- 

To evaluate Pno{t), we use the method of the delay 
function [24]. This yields the probability Pno{t) as the 
square of the norm of the unnormalized state vector: 

Pnoit) = \\\iJno{t))f={lpno{t)\^no{t)) 

= \a{tf + \m\' + \lit)\' + \S{t)\'. (18) 

From Eqs. (flT)) and (fTS]) one obtains 

m = Hnoit)){^no{t)\ + \e{t)\'\e){e\, (19) 

where we have defined 

\e{t)f ^ Py,,{t) ^ 1 ^ Pnoit) . (20) 

The quantities \a{t)\^, |/3(i)P, |7(t)P, \S{t)\^, and |e(t)|2 
represent the probabilities that at time t the system can 
be found either in |a), \b), |c), \d), and |e), respectively. 
To determine a(i), (3{t), -f{t), and (5(t), we have to 
solve the nonunitary Schrodinger equation (jlOp together 
with pT|) . This leads to the inhomogeneous system of 
differential equations, 

f d(t) = -lAXi) - igam , 
/3(i) = -zgaa(i) - (/Ca/2 + ^A,)/3(i) , 
7(t) = -^A',7(i) - Z5b<5(t) , 
t m = -*5ft7(0 - C^b/^ + iAb)d{t) - ^EJTbe^^Pit) . 

(21) 
The differential equations for a{t) and P{t) can be solved 
independently from those for 'y{t) and 6{t). For the initial 
conditions a(0) = 1 and /?(0) = 0, i.e. at time t = the 
atom A is in the state \la) and the cavity A in the vacuum 
state, and defining 



A„ 



{lCa + 2lAa) 



-4g2-^(/C, + 2^A,)A^-A^2 



(22) 
we can write the solutions for a{t) and /3(t), similarly as 
done in [SJ, as 



ait). 



[Kg + 2^A,)/2 - zA;, ^.^^ /iU 
Aa ^'"^ V 2 



cosh 



Kt 



,-[(K;a+2iAa)/4+iA^/2]t 



/?(t) = _^sinh f^^) g-[(C.+2.A„)/4+.A;/2]t 



(23) 



Inserting now in the inhomogeneous pair of differential 
equations for 7(i) and 5{t) the solution obtained for p{t), 
we can determine the solutions for 7(f) and 5{t) with 
the method of the fundamental matrix. For the initial 
conditions 7(0) = and 5(0) = 0, i.e. at time t — Q 
the atom B is in the state |0b) and the cavity B in the 
vacuum state, and defining 



A, 



_ UKb + 2iAb 



z(/C6 + 2zA6)A;-A', 



2 

b ' 



(24) 



we get 

lit) = 9b {f+{t)[g.{t) + h+{t)]-f^{t)[g+{t) + h^{t)]} , 



5{t)=i 



)Cb + 2iAb A' , Ab 



JCb + 2iA 



b_.^_^ 

4 * 2 2 



f-{t)[9+it) + h~{t)] 
f+it)[9-it) + h+it)].i25) 



Here we have defined, for notational convenience, 



„ ,^^ ^ 9aVl^a'^be'' ^[-(K,+2»At,)/4-iA;./2±At,/2]t 
^ ^ A,Ab 



and 



5±(0 



/i±(t) = 



g[(Aa±A6)/2-r-ie]t _ ;^ 
(Aa±A6)/2-T-ze' 

g-[(A„±A6)/2+T+ie]t _ I 

{Aa±Ab)/2 + T + ie ' 



(26) 



(27) 



(28) 



whereT = (/Ca-/Cb+2iAa-2iA6)/4, ande = (A;^-A;)/2. 
In the case of equal parameters for the two subsystems 
A and B, the solutions ([^5]) simplify as 



7(0 



<5(i) 



Kg'e 



2pi0 



A3 



-At^^^_-|^1 g[-(K;+2'iA)/4-iA'/2+A/2]t 



2^id 



ng-e 



A3 
inge 



[eA*_Af-l] 



,[-(K;+2iA)/4-iA72-A/2]t 



i0 



A3 



/C + 2iA .A' A 
4 *T^2" 



[p/* - At - 1] 



X e 



[-(/C+2iA)/4-jA72-A/2]t_ ^'^9^ ' 

A3 



/C+2iA .A' 



[^ 



-At 



-At-1 e 



[~(/C+2iA)/4-iA72+A/2]t 



(29) 



where we have used lim2,^o{[exp(±a;t) — l]/a;} = ±t, 
and defined k = Ka — Kb, K = K-a = ICb, g — (ja — 9b, 

A^Aa = Ab, A' = A;. = A;„ and A = A„ = Ab. 



Using the solutions given by Eqs. ([23)1 and ([251), ^^ 
(|29)) . one can plot the functions |a(i)P, |/3(i)P, |7(0P, 
and |<5(t)P, i.e. the occupation probabilities of the states 
I a), 1 6), |c), and |d), respectively. In Fig. [3] we show 
these probabilities for the case of equal parameters for 
the two subsystems A and B, with g/K. = 10, fl/IC — 10, 
A//C = 1000, and k/K. — 0.9, i.e. the absorption or 
scattering by the cavity mirrors is 10% of the total cavity 
decay. Note that the phase factor e"^ does not play any 
role in the functions considered here. From the figure 
one can see how the dynamical evolution of the source 
subsystem A drives the target subsystem B. Of course, 
for ICt^l, these propabilities are all tending to zero, due 
to the fact that, sooner or later, a photon is absorbed or 
scattered by the cavities mirrors, or is emitted into the 
radiated field, so that the state vector of the system is 
projected onto the state |e). 
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FIG. 3: The probabilities \a{t)\^ (dotted line), |/3(t)P (dot- 
dashed line), |7(i)P (solid line), and |5(i)p (dashed line) are 
shown for the case of equal parameters for the two subsystems 
A and B, where g/K. = 10, fi/ZC = 10, A//C = 1000, k//C = 
0.9, and 2(7 + 7')/A < 1. 



FIG. 4: The concurrence C[psx(t)] between the two atoms is 
shown for the case of equal parameters for the two subsystems 
A and B, where g/K. = 10, ^/IC = 10, A/X: = 1000, 2(7 + 
7')/A < 1, and k/ZC = 1 (dashed line), k/K. = 0.9 (solid 
line), n/K, — 0.8 (dotted line). 



III. UNCONDITIONAL PREPARATION OF 

ENTANGLEMENT 

In this section we study the dynamical generation of 
the entanglemet between the two subsystem A and B. In 
particular the entanglement evolution will be analyzed 
by means of the concurrence. We will also see that the 
entanglement generated between the two atoms can be 
stored by switching off the Raman interaction. 



A. Entanglement evolution 

For the situation under study, the two atoms constitute 
a pair of qubits. An appropriate measure of the entan- 
glement for a two qubits system, often considered in the 
context of quantum information theory, is the concur- 
rence [35| . Given the density matrix p for such a system, 
the concurrence is defined as 

C(p)==max{0, v^- V^- V^- V^} , (30) 

where Ai > A2 > A3 > A4 are the eigenvalues of the ma- 
trix p = p{(Ty ® ay)p*((jy ®(Jy). Hcre Gy IS thc Pauli spin 
matrix and complex conjugation is denoted by an aster- 
isk. The concurrence varies in the range [0, 1], where the 
values and 1 represent separable states and maximally 
entangled states, respectively. 

To derive an expression for the concurrence between 
the two atoms, let us consider the density operator that 
describes the system. It is obtained from the density 
operator /5(i), Eq. (|19p . by tracing over the intracavity 
field states for the two subsystems, Pat{t) — Trcav [p{i)]i 
and is given by 

Pat(t) = |aWni,0)(l,0| + |7(i)nO,l)(0,l| 

+ a(i)7*(t)|l,0)(0, 1| + a*(<)7W|0, 1>(1, 0| 
+ {l-[|a(i)P + l7WP]}|0,0)(0,0|. (31) 



Considering the 4x4 density matrix Pat(Oi related to 
the density operator in Eq. ((3T|) in the two-atom basis 
{|0, 0), |0, 1), |1,0), |1, 1)}, it is easy to show that the con- 
currence C[/3at(0] is, using Eq. ([50)1 . given by 



C[patW]=2|a(i)||7WI- 



(32) 



To analyze the time dependence of this concurrence, let 
us consider the case of equal parameters for the two sub- 
systems A and B. Inserting the analytical solutions ([^5]) 
and ((291) into Eq. ([32]), we show in Fig. g] the function 
C[pa.t{t)] for the parameters g/JC — 10, Vl/JC = 10, and 
A//C = 1000, for different values of k//C. Since the con- 
currence contains only absolute values, the phase factor 
e"^ does not play any role here. From this figure one can 
clearly see that the initially disentangled atoms become 
entangled. In particular, a maximum value for C[pat(0] 
is found for i ~ 28.32//C, where, for the shown cases, 
C[pat(i)] =i 0.73, 0.66, and 0.59. Note that the effects 
due to the absorption or scattering by the cavity mirrors 
are not negligible. For example, the relative variation of 
the concurrence is approximately 10% between the case 
k//C = 1 (no absorption or scattering) and k/K, — 0.9, 
considering the peak at t. Of course, for /Ci ^ 1, the two 
atoms become again disentangled due to the emission 
of the photon in one of the three decay channels. This 
is in agreement with the fact that the release of a pho- 
ton into the environment destroys any entanglement, pro- 
jecting the two-atom subsystem into the separable state 
|0,0). The inclusion of the very rare spontaneous emis- 
sions would only speed up somewhat the decay of the 
entanglement. 

Finally, we note that the concurrence between the two 
intracavity fields can be obtained as well. Let us consider 
the density operator that describes the system of the two 
intracavity fields A and B, obtained from the density 
operator /5(i), cf. Eq. HH), by tracing over the atomic 
states of the two subsystems, Pca.v(t) = Trat [p(0]- It is 



given by 

/5cav(i) = |/3(i)l'|l,0)(l,0| + |5(i)n0,l)(0, 1| 

+ mS*{t)\l,Q){0,l\+ f3*{t)S{t)\0,l){l,0\ 
+ {l-[\m\' + \m\']}\0,0){OA- (33) 

Considering now the 4x4 density matrix 
Pcavit) in the two intracavity- fields Fock basis 
{|0,0),|0,1),|1,0),|1,1)}, the concurrence C[p,,^{t)] is 
given by 



C[peav(i)]= 2 1/3(01 |<5(i)| 



(34) 



Note that the concurrence for the intracavity fields is of 
the same form as the one for the two atoms, cf. Eq. ([5^ . 
when replacing a{t) and 7(i) with f3(t) and S(t). 



B. Storage of entanglement 

We analyze now the possibility to store the entangle- 
ment between the two atoms. As in the previous sub- 
section, let us indicate with i the time when C[pat(i)] 
reaches its maximum value. We consider the case when 
at time t — twe switch off the two lasers, i.e. ^la^i^b — O, 
so that the Raman coupling vanishes. For t > t the in- 
homogeneous system of differential equations becomes 

d(i) = 0, 

7(0 = 0- 

S{t) = -{ICb/2 + ^Ab)<5(t) - ^l^^e'^(3{t) . 

It is immediate to write the solutions for a(t), 7(^)7 and 
P{t) for i > i as 

ait) = ait) , m = /3(i)e-('C"/2+^A„)(t-t) ^ 

lit) - lit) ■ (36) 

Using the solution for /?(<) in the differential equation for 
S{t), one gets, for t >t, 



Sit) 



_ -{ICb/2+iAb){t-i} 



Sit) ~ ^J^JTbe'^Pit) 



,[-(/C„-/C6)/2-i(A„-Ab)](t-t) _ I 
-ilCa-ICb)/2-liAa-Ab) 



(37) 



For equal parameters in the two subsystems Eq. p7l 
reads as 



Sit) 



_ -(K/2+iA){t-i) 



Sit)-Ke''l'(3it)it-t)], (38) 



where we have used lima;^o{exp[a;(t — i)]~l}/a^= (^^f)- 
Using these solutions it is clear that, for t > t, the con- 
currence between the two atoms remains constant and its 
value is 









^-, 







20 .^ 30 

JCt 



FIG. 5; The concurrence C[pa.t(t)] between the two atoms 
(solid line) and the concurrence C[pcav(i)] between the two 
intracavity fields (dashed line) are shown for the case of equal 
parameters for the two subsystems A and B. The parameters 
are g/lC = 10, 0//C = 10, A//C = 1000, and «://C = 0.9. At 
i — 28.32/ K, the two lasers are switched off, so that fi = for 
t> i. 



C[pat(t)] = C[pat(t)] = 2|a(0||7(i)l- 



(39) 



This expression shows that the entanglement between the 
two atoms can be stored even for times with /Ci 3> 1. 
This situation could be realized, for example, by using a 
hyperfine transition in ^Be ions, whose coherence time 
between the two internal levels was reported to be several 
minutes [3g|. The experimental setup could be similar 
to those in Refs. [33, [33|. Note that the concurrence 
between the two cavity fields, given by Eq. ([34)) . is instead 
decreasing, and is quickly vanishing. This is due to the 
fact that a photon in the cavity is, sooner or later, either 
emitted in the radiated field, or absorbed or scattered by 
the cavities mirrors. 



The behavior for the two concurrences is shown in 
Fig. [HI where we have used for the two atom-cavity sub- 
system the same parameters as in Fig. [31 with k/JC = 0.9. 
In this case the two laser beams are turned off at time 
i = 28.32//C, where the concurrence attains its maximum 
value of C[pat(f)] — 0.66. Note that one gets i ~ A/igfl). 
This justifies, in agreement with the discussion in Sec. [Ill 
that one can omit the effects of spontaneous emissions in 
the time interval [0,i]. Moreover, also for t > i sponta- 
neous emissions are negligible. In fact, when the two 
lasers are switched off, the only possibility to have a 
jump related to spontaneous emissions is via the cav- 
ity coupling, i.e. proportional to the terms gaa and gbb. 
This contribution is negligible not only because of the 
large detuning iTa,T'^,Tb,T'ij ^ 1), but it is also van- 
ishing because the cavities are, for t > i, practically in 
the vacuum state. For example, with the values used in 
Fig. [51 already a,t t = t the two decaying functions /3(i) 
and Sit), cf. Eqs. ([M]) and ([HT]) . have negligible values, 

\m?^\sii)\'^om. 



IV. CONDITIONAL PREPARATION OF 

ENTANGLEMENT 

Let us now turn our attention to the case of a condi- 
tional preparation of the entanglement between the two 
atoms, and its subsequent storage. This new situation 
is obtained by introducing a photodetector of quantum 
efficiency rj that monitors the radiated field, as indicated 
in Fig. [TJ We are interested to study the case when "no 
click" occurs, i.e. a conditional evolution under imper- 
fect detection. When a "click" at the photodetector is 
recorded, the conditional preparation is not successful, 
and the preparation procedure has to be repeated again. 

In order to properly treat this problem, let us introduce 
the following consideration, cf. [39] and Appendix „ As 
it has been already mentioned in Sec. |TT1 the probability 
for a jump Ji to occur in the time interval [t, t + dt) is 
given by Pi{t) = {JjJi)tdt. The increment in this time 
interval for pycs{t), cf. Eq. ([20]) . is equal to 



?t 



dpyes(i) = {JlJi)t dt + {JlJ2)t dt + {Jlh)t dt . (40) 

Using Eqs. (O and (O one obtains, by integrating 
Eq. (gOD, that 



Pycs(t) =Prad(i) +Pedys{t) 



(41) 



where 



P..d{t)^f dt'{jlj,)t'^Kaf dt'm')\^+Kj dt'\S{t')\^ 

Jo Jo Jo 

+ 2 4^jr^ f dt'Re[p*{t')S{t')e~''^] , (42) 

Jo 



'/ ft L \,, I / rli'l ft 



and 



Pabs(i) = / dt'{JlJ2)t' + / dt'{JlJ^)t' 

JQ Jo 

= n'['dt'm')\^ + K',fdt'\5{t')\\ (43) 
Jo Jo 

The function prad(i) represents the probability that a 
photon is radiated by the cascaded system in the time 
interval [0,i], and Pabs(i) the probability that a photon 
is absorbed or scattered by the cavity mirrors in the same 
time interval. Note that because 5{t) contains an overall 
factor e*"^, cf. Eqs. pS)) and (|26p. the phase (j) is irrelevant 
in Eq. (|45)l . 

Let us now assume that somehow we know that for 
sure in the time interval [0, t] a photon has been released 
by the cascaded system into its environment. In this case 
one would have that Pyes(t) = 1, and Pno{t) = 0, i.e. we 
know for sure in which of the two possible realizations 
the system is found at time t. The density operator p{t) 
that describes the system is then given, cf. Eqs. (|16p and 
(fT7)l . by p{t) = Pyesit) = |e)(e|. It follows that the two 
atoms are in the separable state |0, 0), and, obviously, 
the related concurrence is equal to zero. The release of 



a photon in the environment destroys any entanglement 
between the two atoms. 

If we now assume that we are in the opposite case, i.e. 
that somehow we know that for sure in the time interval 
[0, t] a photon has not been released by the cascaded sys- 
tem into its environment, then pycs(i) = 0, and Pno(^) = 1- 
The density operator p{t) that describes the system is 
given, in this case, by p{t) = Pno{t), cf. Eq. (fH)). The 
reduced density operator of the system consisting of the 
two atoms, Pat|no(0 = Trcav [Pno(i)], is now given by 

Pat|no(i) = ^--||«(i)ni,0)(l,0| + |7(i)nO,l)(0,l| 
Pno(i) 1^ 

+ a(i)7*(t)|l,0)(0,l|+a*(t)7(t)|0,l)(l,0| 
+ [|/3(t)p + |5(t)p]|0,0)(0,0|}, (44) 

where Pno{t) is given by Eq. ITSl) . It is easy to show that 
from Eq. ([50)1 the concurrence C[pat|no(0] is equal to 



C[pat|„o(i)] 



2\a{t)\\j{t)\ C[pat(t)] 



Pno{t) 



Pno{t) 



(45) 



where we have also used Eq. ([32]) . Because Pno{t) < 1, 
it follows that C[pat|no(i)] > C[pat{t)]- In other words, 
the knowledge that no photon has been released by the 
cascade system increases the entanglement between the 
two atoms. 

Let us now consider the case when a photodetector of 
given efficiency rj is used to monitor the radiated field. 
It is possible to show, by using the quantum trajectory 
method J32, i40|] , see Eq. (|A.11[) in the Appcndixn that the 
probability of not recording a click at the photodetector 
up to time t is given by 

Po{t) = Pno(t) + (1 - '7)Prad(i) +Pabs(i) = 1 - ?7Prad(i) 
= 1 - r; + 77 [pno{t) + Pabs(0] , (46) 

where we have also used Eqs. (|^T|) and pp)) . The con- 
ditional state given that the detector does not record a 
photon, is a weighted sum over the conditional density 
operators reached via the various records of this type, 
i.e. null-measurement at the detector that monitors the 
radiated field. This yields, cf. Eq. (|XT2l) . 

Po{t) = — 7TTlPno(i)/5no(i)+[(l-'7)Prad(0 

+ pMt)]Pycs{t)} = ^TT|lV^no(i)>(^no(t)| 

+ [{l-v)Pr.d{t)+pMt)]\e){e\} , (47) 

where pno{t) and pycs{t) are given by Eq. (|14p and 
Eq. (fT6|) . respectively. Note that if 77 = we obtain, 
from Eq. (|i^ . that po{t) = 1. In this case we return 
to an unconditional evolution, and pait) becomes again 
pit), cf. Eq. (gT]) and Eq. dH]). For r/ = 1, and no photon 
absorption or scattering by the cavity mirrors, i.e. when 
Pahsit) = 0, we obtain, from Eq. (|46p. that po(t) = Pno- 
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This is the case where we can be sure that no photon 
has been lost by the system, and po{t) becomes Pno{t), 
cf. Eq. gll) and Eq. pj]) . 

For the reduced density operator of the system con- 
sisting of the two atoms we have that 

/5at|o(0 = Trcav [/5o(<)] = ^^^Pat|no(i) 



[(l-ry)prad(0+Pabs(0] 



0,0)(0,0|, (48) 



where Pa,t\no{i) is given by Eq. ([44]). Considering now the 
4x4 density matrix Pat\oit), related to the density oper- 
ator of Eq. dlHl) in the basis {|0, 0), |0, 1), |1, 0), |1, 1)}, it 
is easy to show that the concurrence C[/9at|o(^)] is, using 
Eq. jSni), given by 



C[Pat|oW] - 



P0{t) 



C[Pat(i)] 



(49) 



where we have used Eq. (j45p . This is the concurrence 
between the two atoms in the presence of a photode- 
tector, of given efficiency, informing us that no photon 
has been registered in the radiated field. If 77 = 0, we 
have po(^) — I7 and the concurrence is again equal to 
C[pat(0]\ as in Eq. (132D- Note that po{t) < 1, so that the 
concurrence C[pat|o(0] ^ C[pat{t)]. Moreover, po{t) is a 
decreasing function with 77, so that for 77 = 1 it reaches its 
minimum value, and, consequently, cf. Eq. (|49p . the con- 
currence C[pat|o(0] reaches its maximum value. In this 
case, and for perfect mirrors, i.e. for Pabs(i) = 0, one has 
Po{t)—pnoit), and the concurrence is given by Eq. (|l5)) . 

Let us now consider the case analyzed in Fig. [5l but 
with the presence of a detector of efficiency ry that mon- 
itors the radiated field. Because |/3(i)P ^ l^(i)P — 0, 
for /Ci 3> 1, this implies that in this case one has, using 
Eq. (I18l),p„o(t)~|a(i)p-l-|7(i)p. With the values consid- 
ered in Fig. [51 this gives Pno{t) — 0.66. Moreover, from 
Eq. (j43p , and using the same parameters as in Fig. [3 we 
obtain, for ICt ^ 1, the value Pa,hs{t) ~ 0.2 . Consider- 
ing that the single-photon detector-efficiency has alr eady 
reached a value of approximately 77 = 0.88, cf. Ref. [4l| . 
from Eq. (j46p we obtain, for these parameters, the value 
Po{t) ~ 0.88. Because for t > i the concurrence C[pat{t)] 
is given by Eq. ([55)) . i.e. C[pat{t)] ~ 0.66, we obtain from 
Eq. ^ that, for ICt > 1, C[pat|o(i)] - 0.75. This is the 
value of the concurrence stored between the two atoms, 
with an enhancement of approximately 14%. Note that in 
the ideal case of r; = 1 and no photon absorption or scat- 
tering by the cavity mirrors, then po{t) = Pno{t) — 0.66, 
cf. Eq. ([^S)) . so that one would obtain, with the chosen 
parameters, C[pat|o(0] = C'[pat|no(i)] - 1, cf. Eqs. ([45|) 
and ([15)) . This value of the concurrence is related to the 
fact that \a{t)\ ~ \-f{t)\. 

Finally, an important question is related to the prob- 
ability of successfully realizing the whole process of the 
conditional preparation and storage of entanglement. We 
know that if the detector registers a photon, then the 
density operator for the cascaded system is given by 



Eq. ([16)) . and no entanglement is present between the 
two atoms. The probability of a successful realization of 
this scheme is given by the probability that the detec- 
tor does not register any photon, probability given by 
Po(i) = 1 ~ '7Prad(0- If ^ — Oi i-C- whcu we are in the case 
where no detector is present, this procedure is always suc- 
cessful, but the concurrence is not enhanced and remains, 
as in the previous section, given by Eq. ([5^ . Note that, 
because we are interested here in a null-measurement 
conditional preparation, the above considerations remain 
valid also when a click at the photodetector is coming 
from a possible dark count. This only increases the prob- 
ability that the whole procedure has to be repeated from 
the beginning. Using the same parameters as in Fig. [Sj 
and a detector efficiency ry = 0.88, cf. Ref. [41|, we have 
that po{t) ~ 0.88. This means that for approximately 
88% of the cases, the conditional preparation and stor- 
age of the entanglement between the two atoms, with a 
value C[/9at|o(0] — 0.75, is successfully realized. For the 
remaining cases we have to repeat the whole procedure 
from the beginning. 



V. CONCLUSIONS 



The dynamics of a cascaded system that consists of two 
atom-cavity subsystems has been analyzed. Considering 
the two atom-cavity subsystems driven by a Raman inter- 
action, the evolution of the open quantum system under 
study has been described by means of a master equation. 
By using the quantum trajectory method, analytical so- 
lutions for the dynamics of the system have been ob- 
tained. The entanglement evolution between two stable 
ground states for the two atoms, constituting a two-qubit 
system, has been studied using the concurrence. A sim- 
ilar analysis has been performed for the two intracavity 
fields. 

The dynamical evolution of the system shows that the 
two initially disentangled qubits reach states of signifi- 
cant entanglement. Moreover, it has been shown that 
the entanglement generated between the two atoms can 
be stored by switching off the Raman coupling. Subse- 
quently, we have analyzed how the entanglement between 
the two atoms can be enhanced, by monitoring the radi- 
ated field with a photodetector of given efficiency, via a 
null-measurement conditional preparation. 
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APPENDIX 

Let us analyze the problem of the conditional evolu- 
tion under imperfect detection [33, Ho] by considering a 
system with two output channels. Channel 1 is moni- 
tored by a detector of efficiency rj and associated with 
a jump operator J7i. Channel 2 is monitored by a de- 
tector of unit efficiency, which could represent the envi- 
ronment into which the system releases a photon, asso- 
ciated with a jump operator ^2- One can treat non-unit 
detection efficiency by introducing a beam splitter into 
channel f , with transmittivity y/rj. Now the system has 
three output channels, the transmitted and the reflected 
parts fx and fpt, respectively, due to the beam splitter, 
and the original channel 2. In principle, all three chan- 
nels could be thought of being monitored by detectors 
of unit efficiency. In the following we will indicate the 
three channels with i = 1t,1r,2. We are interested in 
the conditional evolution under null-measurement at the 
photodetector that monitors the transmitted beam, i.e. 
channel It. The master equation for the density oper- 
ator that describes the system can be formally written, 
cf. [3l|, as dp{t)/dt — Cp{t), where the superoperator 
L = Lb + S\^ -t- S\^ -t- ^2 is in the usual Lindblad form. 
The between-jump superoperator Lb is given by 



Lb- = -z 



Hs, 



2 Z^ y^i "^^ 



U j \Ji 



(A.l) 



where Hs is the system Hamiltonian, with the jump oper- 
ators for the transmitted and the reflected channels given 
by 



The three jump superoperators are defined as 
5,-=i, -i; (z = 1t,1r,2). 



(A.2) 



(A.3) 



Let the initial density operator be /5(0). Because we 
are interested in a system where we can have at most 
one jump, there are four records of interest. First, the 
record where neither detector clicks in the interval [0,i); 
for it we have the probability and conditional density 
operator |31 | 



^Cst 



P„o(0 = Tr[e^«*p(0)], p„o(t) = 



m 



Tr[e^i5*p(0)] 



(A.4) 



We have then the record given by a photon detected in 
the time interval [t',t' + dt), with t' < t, at one of the 
detectors i {i — It, 1r, 2) and the other two detectors not 
clicking. For it we have the probability and conditional 
density operator (z = It, 1r, 2) 



Pi{t) = Tr[e^«(*-*')5ie^«*'/5(0) dt , 



(A.5) 



The probability for no click at detector It up to time t 
is given by a sum over all events with no click in channel 
It, 

P0{t) = Pno{t) + f dt'pi^it') + [ dt'p2{t') . (A.6) 
Jo Jo 

The corresponding conditional density operator is a 
weighted sum over conditional density operators of the 
form 



Po{t) 



1 



Po{t) 



Pnoit)pnoit)+ / dt'pi„{t')pi„{t') 



+ / dt'p2{t')p2{t') 



(A.7) 



Following the quantum trajectory method [31], when no 
jump occurs, the system evolves between time io = and 
t via 



pUt) = e^-'m 



(A., 



where Pno(^) i^, i'^ general, not normalized. The evolution 
governed by Eq. (jA.8[) is randomly interrupted by jumps. 
If a jump occurres at time ij, ij G [0,i), the density 
operator collapses into 



S^pLit.^) = J^P'no{ti)Jl ^ |0)(0| [l - 1t,1r,2). 

(A.9) 
Let us now define 



Pi{t) = j dt\jlj^,)r + f dt\jlj^^)t, , 
Jo Jo 



P2{t) = dt'{J^j2)t', 

Jo 



(A.IO) 



where (. . .)t' = Ti[p'^^{t') . . .] . The function Pi(i) rep- 
resents the probability that in the time interval [0, t) a 
photon is emitted by the system into channel It or 1r, 
and P2 (t) is the probability that a photon is emitted into 
channel 2. The probability that the detector It does not 
click up to time t is obtained from Eq. (|A.6[) as 

Poit) = Pnoit) + {l-v)Plit)+P2it) , (A.ll) 



with Pno{t) given in Eq. (jA.4p . Finally, the conditional 
density operator given that the detector It does not click 
is given, using Eq. (|A.7[) . by 

Poit) = ^{pLit)+[{l'V)Piit) + P2mO){0\}.iA.l2) 
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